Renormalization group of probabilistic cellular automata with one absorbing state 
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We apply a recently proposed dynamically driven renormalization group scheme to probabilistic 
cellular automata having one absorbing state. We have found just one unstable hxed point with one 
relevant direction. In the limit of small transition probability one of the cellular automata reduces 
to the contact process revealing that the cellular automata are in the same universality class as that 
process, as expected. Better numerical results are obtained as the approximations for the stationary 
distribution are improved. 



I. INTRODUCTION 



Recently, a general dynamical renormalization group 
(RG) scheme p] has been proposed in order to treat 
nonequilibrium critical phenomena. The method, called 
Dynamically Driven Renormalization Group (DDRG), 
has been applied to self-organized critical phenomena, 
specifically to critical height sandpile models Q and for- 
est fire models §] . Its phenomenological approach takes 
into account the nature of self-organized systems through 
an attractive fixed point. The scheme also provides nu- 
merical values for the critical exponents which are close 
to the ones obtained using computer simulations. 

The scheme consists in coupling a real space RG 
scheme to a stationary condition that drives the RG 
group equations through the parameter space. The sta- 
tionary conditions, involving the stationary distribution, 
must always be approximated. While the authors used a 
stationary probability distribution that neglects any cor- 
relation among different sites, they have mentioned |l]] 
the fact that using more refined approximations should 
improve the values of the critical exponents. 

While the DDRG scheme is quite general, self- 
organized critical systems are special since they have a 
well defined time scale separation (dissipation events be- 
ing instantaneous with respect to the driving) . This pre- 
vents proliferation effects in the real space RG, making 
the calculation of critical exponents easier. 

In this work, we implemented the DDRG scheme to 
another class of systems namely the class of driven diffu- 
sive systems with an absorbing state (the vacuum state) . 
Due to the presence of one absorbing state this class of 
systems is in the same universality as the directed perco- 
lation model H ||. More precisely we have considered 
one-dimensional probabilistic cellular automata with one 
absorbing state. We have treated two types of models. 
One of them is a two-state probabilistic cellular automa- 
ton (that can be interpreted as directed percolation in 
two dimensions) that includes, in the limit of small tran- 
sition probabilities, as particular cases, continuous-time 
processes with one absorbing state such as the contact 
model H and others 0. The other one is a four-state 
probabilistic cellular automaton that includes the model 



introduced by Grassberger and de la Torre g as a par- 
ticular case. 

The appropriate real space RG parameter space for 
these nonequilibrium models is the space of the transi- 
tion probabilities, instead of being the space of coupling 
constants as is the case of equilibrium models, defined by 
a Hamiltonian. In the case of the models studied here, 
with one absorbing state, the RG should be appropriate 
to preserve the vacuum state along the RG trajectory in 
this parameter space. By using a block renormalization 
to treat properly the absorbing state we have figured the 
value of the critical exponent corresponding to the diver- 
gence of the spatial correlation, u± , using three different 
approximations which consider correlations among clus- 
ters up to 1, 3, and 5 neighboring sites in the lattice, 
respectively. Our best calculations give v±_ = 1.04 ± 0.02 
which is rather close to the one obtained from numerical 
simulations reported in || namely = 1.067 ± 0.005. 



II. TWO-STATE MODEL 

The first model is a one-dimensional cellular automa- 
ton with just two state per site. Each site can be either 
vacant, <Ji — 0, or occupied by a particle, Oi = 1. At 
each time step the state of a certain site will depend 
only on the previous states of that same site and its 
nearest neighbors. We consider the most general tran- 
sition probabilities that are homogeneous and symmetric 
in space. The transition probability M^(cr|cr') from state 
a' = (a[, a' 2 , ■ ■ ■ , <j' n ) to state a = (ai, • • • , on) will 
be given by the product 



N 



W(*\<T') = l[w(<T i \*' i _ 1 ,<T , i ,<T , i+1 ), 



(1) 



where N is the number of sites and w(oi\o' i _ 1 ,<j' i ,<j' i+1 ) 
is the one-site transition probability given by the rules 
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The rule w(0|000) = 1 implies that the vacuum state 
is indeed an absorbing state. When P5 = Pi = pz we say 
that the annihilation of particles is spontaneous. Sup- 
pose moreover that pi — 2p\ and that the parameters pi 
and £»3 are very small. In this case the system remains 
most of the time in its previous state. We expect, there- 
fore, that the properties of the present two-state model, 
in the limit p\ — > and p% — > 0, with the ratio p\jpz = A 
fixed, be identical to the contact process with a catalytic 
transition rate equal to A and a transition rate for sponta- 
neous annihilation equal to unity. If p 2 = p\ , and taking 
the same limit, the properties will be identical to a model 
introduced by Dickman (Model A) [| . 

The model is represented by a set of five parameters, 
Pi, Pi, pa, Pi, and ps, which constitutes, as we shall see, 
the RG parameter space. The RG scheme will be con- 
structed in a way that the RG trajectory will be confined 
to this space. 

III. FOUR-STATE MODEL 

Grassberger and the la Torre model M is defined as 
follows. Each site of a one-dimensional lattice is either 
occupied by one particle or it is void. At a certain time 
step the state of the system may be defined by the vec- 
tor er = (o"i, 02, (T3, • • • , <tn) where <Ji = or 1 according 
whether the site i is vacant or occupied. In each time 
step, all sites are updated in two stages. 

1) In the first stage each particle is spontaneously an- 
nihilated with probability c. 

2) In the second stage every surviving particle will 
generate, with probability p, a new (unique) parti- 
cle, which will be placed in one of its nearest neigh- 
boring sites, randomly chosen. In other words, for 
each site with a particle, we chose a neighboring 
site with probability p/2. If the site was originally 
void it becomes occupied, and if it was occupied 
it remains as such. We have modified slightly the 
original model by introducing the parameter p. The 
original model of Grassberger and de la Torre is re- 
cast when p = 1. 

Defined in this way, the transition probability Waaler') 
from a state er', to another state <r, can not be written 
as a product of independent transition probabilities as- 
sociated to each site w(cri|(T^_ 1 , a[, cr' i+1 ), as in ordinary 
cellular automaton. However, if we enlarge the number 
of states in each site by introducing three types of parti- 
cle, then it is possible to map the model into a four-state 
cellular automata. This mapping is outlined in Appendix 
A. 

The four-state probabilistic cellular automaton equiv- 
alent to Grassberger and de la Torre model is defined as 
follows. Each site of a one-dimensional lattice can be ei- 
ther empty (E), <Ji = 0, or occupied by a neutral particle 
(N), crj = 1, or by a rightist particle (R), Oi — 2, or by a 



leftist particle (L), cr, = 3. At each time step, every site 
of the lattice is independently updated according to the 
rules. 

1) If the site is occupied by one particle of any type 
N, R, or L, then one out of four possible events 
will take place. 

a) The particle is annihilated, that is, the site 
becomes empty, with probability c, or 

b) becomes a particle of type N with probability 
a, or 

c) becomes a particle of type R with probability 
6/2, or 

d) becomes a particle of type L with probability 
6/2. 

Here a = (1 — c)(l — p) and 6 = (1 — c)p. 

2) In case the site is empty (state E) one has to look 
to its neighboring sites. 

a) If its left neighbor is type R or its right neigh- 
bor is of type L, the site remains as E with 
probability c, becomes N with probability a, 
becomes R with probability 6/2, or becomes 
L with probability 6/2. 

b) If on the contrary, its left nearest neighbor is 
not a particle of type R and is right nearest 
neighbor is not of type L, the site remains va- 
cant. This last rule implies that the state with 
all sites empty is an absorbing state. 

Transition probability W(<r|<7') from state a' = 
(a[, a' 2 , • • • , <j' n ) to state a = (<ri, 02, • • • , <jn) can be writ- 
ten as the product 

N 

W(a\a') = l[w(<T i \a , i _ 1 ,<T'i,<Ti + i), (3) 

i=l 

where N is the number of sites of the lattice. The one- 
site transition probability w(ai\a' i _ 1 ,a' i ,a' i+1 ) is written 
down in the Appendix A. 

In order to apply the RG scheme we enlarge the space 
of parameters but preserve the existence of the absorbing 
state. Demanding besides that the rules should be homo- 
geneous in space and invariant by exchanging the states 
R and L, the transition probability w(ai\a' i _ 1 , a[, 
will be defined in the most general way by 59 parame- 
ters. We call such a probabilistic cellular automaton the 
four-state model. 
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IV. RENORMALIZATION SCHEME 



and equations (||) , (0) , and @ . We obtain [|J 



Here we use a real space RG scheme [p|-[l 1[ which renor- 
malizes the transition probability W. The succession of 
RG transformations corresponds to a trajectory in the 
space spanned by the parameters that defines W. The 
scheme we use is an implementation of the DDRG 
and is accomplished by transforming cells of b sites into a 
cell of just one site. To treat the vacuum state properly 
any cell with at least one particle renormalizes into an 
occupied site. Only cells with no particles renormalizes 
into a vacant site. 

Let 1Z(t\(j) be a condition probability of state r given 
the state a with the following properties 



n(r\a) > 0, 



i. 



(4) 



The vector a = (<ji, 02, • • • , crjsr) represents the state 
of a system with N degrees of freedom and the vector 
t = (t\ , T2, • • • , tjv' ) represents the state of the renormal- 
ized system with N' = N/b degrees of freedom, where b 
is the size of the renormalization block. 

Let .P n (<7, a') be the probability of occurrence of state 
a 1 at a given time and state a at n time steps later, that 
is 



P n (a,a') = W n (a\a')P(a') 



(5) 



where P(cr) is the stationary probability distribution 
which satisfies the equation 



P(a)=J2w n (a\a')P(a') 



(6) 



for any value of n, where W n (a\a') is the transition prob- 
ability from state a' to state a in n time steps. Similarly, 
for the renormalized system, let P(t, t') be the probabil- 
ity of occurrence of state r' at a given time and state r at 
one time step later. The RG transformation are obtained 
by demanding that || 

p (r, r') = E E nrW)n(r' \a')P n (a, a') (7) 

a a' 

from which follows 



since 



(8) 



(9) 



To get the desired renormalized transition probability 
W(t\t') we use 



W(t\t') 



P(t>) 



(10) 



This equation, however, is not properly a transformation 
between the transition probabilities W and W, since the 
yet unknown stationary probability P(cr) appears in the 
right hand side of equation (11). However, if we use the 
balance equation (^|) for P(a) then a closure condition for 
the renormalization group is obtained. The closure con- 
dition plays the role of the driving condition that forces 
the system to be in the stationary state at each step of 
the transformation. In this sense the present DDRG may 
be thought of as a renormalization of the stationary state. 

At each state of the transformation the transition prob- 
ability W always describes an irreversible process so 
that the corresponding stationary solution P(a) obtained 
from the closure condition ([|) will not be related, a priori, 
to a Hamiltonian, that would be the case if the process 
obeyed detailed balance. In this way the present RG 
scheme is distinct from the ordinary real space RG used 
in equilibrium systems in which the parameters of the 
Hamiltonian are renormalized. 

The closure relation (^|), however, cannot actually be 
solved so that approximations should be used. Here we 
have used three different approximations which consider 
correlations among clusters up to 1, 3, and 5 neighboring 
sites. Equation ( |ll| ) together with a given approximation 
provides then a well defined RG transformation W — * W . 

Assuming that the renormalized transition probability 
can also be written as a product of independent transi- 
tion probabilities, that is, 



N' 



W(r\r') = l[w(r k \4_ 1 ,r'kA + x), 



(12) 



fc=i 



one can write the RG equations that relates the old one- 
site transition probability w(ai\a' i _ 1 , a[, to the new 
one wiTklT^T^T'^). 

We have used only renormalization cells with size b = 2 
and chosen 1Z in the form 



N/2 



Tt{T\<j) = R(Tk\(T2k-l,&2k), 



(13) 



fc=l 



with 



R{T k \a 2 k-i,o-2k) > 0,and ^ R( T k\<?2k-i, c 2 fc) = 1 

(14) 

To preserve the absorbing nature of the vacuum state we 
have chosen R with the properties 



i?(0|0,0) = 1 



(15) 
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and 



R(0\a 2k -i,<T2k) = 



(16) 



whenever er 2 fc-i ^ 0, or a 2 k ^ 0. 

For the four-site model we used the following values 

-R(OjOO) = 1 i?(l|01) = l i2(l|10) = 1 J2(l|ll) = 1 

i?(2|02) = 1 i?(l|03) = 1 i?(l|20) = 1 R(3\30) = 1 

fl(2|12) = l i?(l|13) = l i?(l|21) = 1 R(3\3l) = 1 

J2(2|32) = 1/2 E(3|32) = l/2 i?(2|22) = 1 i?(l|23) = 1 

The other matrix elements of R are zero. The first ma- 
trix element has been chosen in order to preserve the 
absorbing nature of the vacuum state and the rest of the 
elements where assigned bearing in mind the physical pic- 
ture that an R particle will give rise to a new particle to 
the right in the next time step, while an L particle will 
generate a new particle to the left (N particles will not 
generate new particles). 



V. RENORMALIZATION ALGORITHM 

The temporal coarse graining will be done using two 
time steps that is n — 2. Using the equations of the pre- 
vious section we can write down the equation that relates 
w to w for the case n = 2 (see figure [l]) . 
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FIG. 1. Diagram showing the blocking scheme procedure. 
Numbers correspond to the indexes used in equation ([n]). 

It is given by 

w{t 4 \t 2 ,t 3 ,t 4 ) = R(n\aia 2 ) X 

<7\<J<2<77 ■■■<7\<2 

T{a\,a 2 \a 7 , ■ ■ -,a 12 )p(a 7 , ■ ■ ■ ,a 12 \T 2 ,T 3 ,r 4 ), (17) 

where 

T(a 4 a 2 \a 7 , as, erg, a w , an, a 42 ) — 

'^2w{ai\a 3 ,a 4 ,a 5 )w{a 2 \a 4 ,a 5 ,a 6 ) x 



&3 <7 4< 7 5 a '6 



w(a 3 \a 7 ,a 8 ,a 9 )w(a 4 \a 8 ,a 9 ,a w ) x 
w(a 5 \as,a w ,an)w(a & \a w ,an,Oi 2 ) , (18) 



p(a 7 ,as,a 9 ,a w ,an,a 12 \T 2 ,T 3 ,T 4 ) = [P(r 2 , t 3 , t 4 )] x 
R(T 2 \a 7 ,as)R(T 3 \a 9 ,a w )R(T 4 \an,^i 2 ) x 
P(&7, &8, 09, cio, an, C12) , (19) 

where 

P(t 2 ,t 3 ,t 4 ) = s ^R{T 2 \a 7 ,a 8 )R(T 3 \a 9 ,a w ) x 

R(T 4 \an,a 12 )P(a 7 , a 8 ,a 9 ,a w ,an,<Ji 2 )- (20) 

Here the subscripts refer to the site numbers appearing 
in figure [l]. 

In order to solve the system of equations (jl7j-p0"|) one 
must resort on approximate methods to estimate the sta- 
tionary weights P(a 7 , • • • , a 42 ). The simplest approxima- 
tion, sometimes known as simple mean field approxima- 
tion, consists in neglecting correlations among different 
sites 



12 



,012) 



(21) 



i=7 



where P{ai) is the solution of 



and 



J2 w^i\^a 4 )P(a 2 )P{a 3 )P(a 4 ). (22) 

However, correlations are actually taken into account 
in the geometrical aspects of the blocking procedure. In 
this way, one obtains non classical critical exponents. 
Given a blocking prescription, the value of the critical 
exponents should improve as correlations are taken into 
account into the stationary probability distribution. In 
order to verify how important the changes will be, we 
have used three different approximations for the station- 
ary distribution. The first one being ( fjl|f22| ), while the 
other two are mean field approximations which consider 
correlations up to clusters of three sites and five sites 
respectively jl2j. 

Due to the number of terms involved in equations ( [l7| - 
|20| ) it is not possible to determine analytically the fixed 
point of the transformation. So, we have performed it 
numerically, using initial values for the transition proba- 
bilities corresponding to the model of interest. In the case 
of the two-state model we start with w(ai\a' i _ 1 , a[, a' i+l ) 
given by (Q) with p§ = p 4 = p 3 and p 2 = 2p 4 . In the case 
of the four-state model we start with w{ai\a' i _ 1 , a'i, a' i+1 ) 
corresponding to the Grassberger and de la Torre model, 
with p = 1. 

In each iteration of the RG, given the set of pa- 
rameters w, one has to find the stationary solution for 
P(a 7 , ■ ■ ■ , C12). This has been done by iterating the time 
evolution equation for the model (using one of the three 
approximations) until reaching convergence. For the one- 
site approximation, given by (|22]), 10 4 iterations were 
enough. As approximations are refined equations become 
highly nonlinear, and for the five-site approximation least 
10 5 iterations are needed. 
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VI. RESULTS 

For the four-state model, the RG equations behaved in 
the following way. For small values of c, the set of tran- 
sition matrix elements flows towards an attractive fixed 
point characterized by c = 0, and a lattice full of parti- 
cles. On the contrary, for values of c high enough, the 
RG equations are driven to a different attractive fixed 
point, this time characterized by c = 1 and a lattice 
without particles. In this case we used only the one-site 
approximation and found a critical value of c given by 
c cr = 0.3568. Starting around this values the representa- 
tive point of the parameter set spends a long time near 
an unstable point before it leaves towards any of the two 
attractive fixed points. Figure 2 shows a projection of 
two trajectories in the parameter space in terms of two 
of these parameters: w(0|010) and I — ii>(0|101). 



1.0 <?c 
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x 

FIG. 2. Projection of two trajectories in the parameter 
space of the four-state model in terms of two of these param- 
eters: x = 1 — tu(0|101) and y — w(0|010). Each trajectory 
reaches a different fixed point. The inset shows an enlarge- 
ment around the nontrivial fixed point. 

In this way, it is found only one relevant parameter. 
Since we are dealing only with stationary properties of 
the model it is reasonable to assume that this parameter 
is associated to the divergence of the spatial correlation 
length and not to the temporal correlation length. So, 
figuring the eigenvalue A associated to that parameter 
we get z/j_ = In 2/ In A. The value measured numerically 
is u ± =0.93 ± 0.005. 

From computer simulations results reported in || one 
can obtain i/± = 1.067 ± 0.005, and the critical value 
of c is c cr = 0.279. The discrepancy between the two 
results are mainly due to the poorness of the one-site ap- 
proximation. By increasing the order of approximation 
the results gets better as we shall see in the case of the 
two-site model. 

We have corroborated that, as one would expect, the 
value of v± for the simpler two-state model and one-site 



approximation is identical as the one previously found. 
Since for refined approximations of three and five sites, 
numerical computations become too imposing, we have 
used these two better approximations only for the two- 
state model. The value found using the three-site approx- 
imation was v± = 0.98 ± 0.01 whereas the value found in 
the five-site approximation was u± — 1.04 ± 0.02, which 
is rather close to the one calculated from numerical sim- 
ulations. Bellow we show the coordinates of the unstable 
fixed point for the three approximations, as well as the 
corresponding value of A 

appr. l pi I p 2 | p 3 [ Pi \ Pb j A 
1 0.22794 0.40381 0.26780 0.22116 0.18245 2.105 
3 0.12107 0.22750 0.14761 0.07870 0.03140 2.025 
5 | 0.07107 1 0.16233 1 0.07330 1 0.04783 1 0.01933 |T950 

Another advantage in using better approximations is the 
improvements one obtains in the critical values of non- 
universal parameters. We have consider the following 
initial conditions for the two-state model: p$ = p^ = P3 
and P2 = 2pi. By varying p\ and p^ we have obtained 
the critical values A c of the ratio A = p\ jp^ as shown in 
the table 
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0.1 


0.01 


0.001 


X c (1-site app.) 


2.8027 


2.0833 


2.0088 


2.0008 


A c (3-site app.) 


3.2365 


3.1205 


3.1703 


3.1757 


A c (5-site app.) 


3.2555 


3.1645 


3.2206 


3.2300 



(23) 



which, in the limit p\ ~ > 0, should be compared with the 
critical value of the contact process A c = 3.299 || . The 
fact that the column associated to p\ = 1 also seems to 
converge to this value is a mere coincidence. 

Reasoning along this line, one may wonder why not to 
use the stationary distribution directly from a computer 
simulation of the model. Making large statistics over con- 
figurations in the stationary regime, one should be able 
to accurately estimate the probabilities of the clusters ap- 
pearing on the left hand side of equations (|^) and (pCj). 
While the idea is in principle right, one can not overcome 
in practice the huge amount of time needed to obtain val- 
ues that are accurate enough. Slight fluctuations in the 
estimated values will make trajectories randomly shift 
their destination towards one of the two attractive fixed 
points, depending on the random seed used in the simu- 
lation. This effect takes place even before any trajectory 
is able to reach the linearized domain of the transfor- 
mation around the unstable point. A way to decrease 
fluctuations is increasing the size of the lattice in which 
one performs the simulation. But the precision one gains 
does not grow faster than ~ 1/yN, where N is the size 
of the lattice. 



VII. CONCLUSION 

We have applied a real space renormalization group 
scheme to a class of driven diffusive probabilistic cellular 
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automata having one absorbing state. Two models hav- 
ing been considered. One of them is a two-state model 
that reduces to the contact process in the limit of small 
transition probabilities. The two-state model can also be 
interpreted as a generic directed percolation in two di- 
mensions. The other is a four-state models that includes 
the model introduced by Grassberger and de la Torre in 
a study related to the contact process. We have found, in 
the RG space of parameters, just one non-trivial unstable 
fixed point with one relevant direction. The existence of 
just this unique fixed point reveals that the probabilis- 
tic cellular automata with one absorbing state belongs 
to the same universality class as the directed percolation 
and the contact process, as expected. 

The implementation of the RG scheme nedeed the cal- 
culation of the stationary probability distribution which 
were performed in several levels of approximations. In- 
creasing the number of cluster size used in the approxi- 
mation improved results were obtained not only for the 
critical exponent v± as well as for the non-universal crit- 
ical quantities. 



VIII. APPENDIX A 

The Grassberger and de la Torre model is a stochas- 
tic process in which at each time step particles are being 
created and annihilated. We imagine the process as a se- 
quence os states A, A', A", A'" , ■ ■ ■ each one being given 
by a vector r\ = (ryi, rj 2 , ■ ■ ■ , ?7at) where rji — or 1 accord- 
ing whether site % is vacant or occupied by a particle. We 
may think of each transition, say A — > A' , as composed 
of three stages with two intermediate states B and C, to 
be defined shortly, so that the whole stochastic process 
corresponds to a sequence A, B, C, A' , B', C , A", B", 
C" , A'" , • • • We will then write the transition probability 
W(A'|A) from state A to A' as given by 

W(A'\A) = W 3 (A'\C)W 2 (C\B)W 1 (B\A). (24) 

B C 

where W 3 (A'\C), W 2 (C\B), and Wi(B\A) are the in- 
termediate transition probabilities related to the three 
stages. 

First stage. In the first stage (A — ► B) of the Grass- 
berger and de la Torre model, each particle is annihilated 
with probability c, so that the probability Wi(.B| A) of the 
transition from A to B is given by 

N 

Wx{B\A) = WW = HMMVi), (25) 



where 



«;i(0|l) = c, wi(l|l) = l-c, 
u;i(l|0) = 0, u>i(0|0) = l. 



(26) 



Second stage. In this stage (B — > C) every particle de- 
cides whether it will generate new particles either to the 



left or to the right in the next step. Each occupied site 
will be labeled according to its decision. A particle that 
decides not to generate particles will be labeled by the 
number 1. A particle that decides to generate another 
one to the right, will be labeled by the number 2 and 
a particle that decides to generate a new particle to the 
left, will be labeled by the number 3. A state of type C 
is then defined by the vector u = (01, a 2 , • • • , cfn) where 
<Tj = 0, 1, 2, 3, so that the transition probability W 2 (C\B) 
of the transition from B to C is given by 



N 



W 2 (C\B) = W 2 (a\r 1 ") = Y[w 2 (a t \^), 



(27) 



i=i 



where 



w 2 (0\0) = 1, 
w 2 (2|0)=0, 
uto(0|l)=0, 
^(2|1) = |, 



w 2 (l\0)=0, 
w 2 (3\0) = 0. 
w 2 (l\l) = l-p, 
«*(3|1) = §, 



(28) 



Third stage. In this state (C — ► A r ), particles are ef- 
fectively created. Each occupied site remains occupied. 
Each vacant site becomes occupied if the site at right 
(left) is occupied by a particle of type 3 (2). 

A configuration of type A' is expressed again in terms 
of the two state variables, rji = or 1, and the transition 
probability W2(-A'|C) of the transition from C to A' is 
given by 

N 

W 3 (A'\C) - W 3 (r?V) = l[w 3 (r 1 l\a i _ 1 ,a i ,a i+1 ), (29) 

where the transition probability w 3 (l\<7i-i, <Ji, <Ji+i) to 
the state 77- = 1 is given by 

w 3 {l\ui-i,ai,a i+ i) = 1, (30) 

if <7j ^ 0, for any value of <7j_i and Uj+i, and 

w 3 {l\o-i-i,0,a i+1 ) = 1 (31) 

if<7j_i = 2 or (Tj+i = 3. In other cases w 3 (l\o~i-i, o~i, Oi+i) 
vanishes. The transition probability w 3 (0\<Ji-i, di, <7,+i) 
to the state ry- = is just given by 

w 3 (0\ai-i,ai,ai + i) = 1 - 103(1(^—1, a,, <7j+i) (32) 

It ease to check that W(A'\A) — W(r)'\rj) cannot be 
written as product of independent transition probability 
of each site as in an ordinary cellular automaton. How- 
ever, the transition probability W C (C'\C) = W c (a'\a) 
from state C to state C can. Indeed, from 

W C (C'\C) = J2J2W 2 (C , \B , )W 1 (B'\A')W 3 (A , \C) 



B> A> 



(33) 



G 



that is from 



And finally, when a\ — 0, and <r' i _ 1 7^ 2 and cr' i+1 7^ 3, 



W c {a'\o)=Y.Y. W ^'W') W ^'W)Wz{v'W) (34) 



77" 77' 



we get 



A' 



^3(^1^-1,^,^+1) 



which can be written in the form 

N 

W c (a'\a) = ]^ii;(CT,-|CT i _i,er i ,erj +1 ) 



(35) 



(36) 



i=i 



where 



w(a[\a l -i,a l ,a l+ i) = ^ E ^(o-^-Owi (v"\Vi) x 

^(lil^-b^^i+i) (37) 

The Grassberger and de la Torre process can then be 
viewed as sequence of states C, C, C", C 1 ",---, each 
one being given by a vector a = (<7i, 02, • • • , <tn) where 
<Ji — 0, 1, 2, or 3 according whether site i is either vacant, 
or occupied by a particle that does not generate another 
particle (neutral particle) , or occupied by a particle that 
generates another one to the right (a rightist particle) 
or occupied by a particle that generates another one to 
the left (a leftist particle). Therefore, it is equivalent 
to an ordinary four-state cellular automaton whose rules 
are defined by equations @, (|f), (f|), and (^tt§. 
From these equations we may write down the transition 
probability w{a' i \ai-i, <Ji, cn+i) in the form 



M0K-i>^>^+i) = c ! 
w(2\a' i _ 1 ,°'i,°' i+ i) = b/2, 



(38) 



i) = V2, 



if o"^ 7^ 0, independently of the states taken by and 
. For the case where a[ — 0, and either cr' i _ 1 = 2 or 
= 3, one has 



'i+i 

T i+1 



w(0\a^ 1 ,0,a' l+1 ) = c, 
w(l|^_ 1; 0,^ +1 ) = a, 
W (2|^_ 1 ,0,a,: +1 ) = 6/2, 
w{2,\a' i _ 1 ,Q,a' i+l ) = b/2, 



(39) 



W (0|a^ 1 ,0,a J ' +1 ) = l, 



(40) 



(1 



The parameters a and b are related to p by a 
p)(l — c) and b = p(l — c). 

The rule w(0|0, 0, 0) = 1 implies that the state with all 
sites vacant is indeed an absorbing state. 
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